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^-^' Abstract 

1-5 ■ A study of ray trajectories was undertaken for the Tamm medium which represents the spacetime 

of a cosmic spinning string, under the geometric-optics approximation. Our numerical studies revealed 
that: (i) rays never cross the string's boundary; (ii) the Tamm medium supports evanescent waves in 
regions of phase space that correspond to those regions of the string's spacetime which could support 
^ ' closed timelike curves; and (iii) a spinning string can be slightly visible while a non-spinning string is 

CO , almost perfectly invisible. 

^ — I , 

T-H . 
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O ! 1 Introduction 

Parallels have been established between exotic optical phenomenons associated with certain metamaterials 
and those associated with curved spacetimcs [1, 2, 3]. For example, negative-phase-velocity propagation 
^^ ' of light — a notable property supported by certain negatively refracting metamaterials [4] — is associated 

with various spacetime metrics in a noncovariant formalism [5, 6, 7, 8, 9]. This parallelism is perhaps not 
surprising given that light propagation in vacuum subjected to a gravitational field is formally equivalent 
to light propagation in a nonhomogcneous bianisotropic medium — called a Tamm medium [10] — in flat 
spacetime [11, 12, 13]. The practical realization of Tamm mediums is edging ever closer due to rapid advances 
in the science of nanostructured materials, especially those relating to metamaterials supporting negative 
refraction or cloaking applications [14, 15]. 

The constitutive relations for the Tamm medium representing a spinning cosmic string were recently 
established [3] . Cosmic strings are topological defects associated with phase transitions which may represent 
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traces of the very early universe, as is comprehensively described elsewhere [16, 17]. In the following sections 
we report on the trajectories of light rays in the Tamm medium representing the spacetime metric of a 
cosmic spinning string, in the geometric-optics regime. In particular, we report on a cloaking phenomenon — 
analogous to that associated with certain metamaterials [18] — which arises when the cosmic string's angular 
momentum is set to zero, thereby rendering it almost perfectly invisible to optical probes. 



2 Quasi— planewave analysis 



We consider a cosmic spinning string that is aligned parallel to the z axis. Let p = \Jx'^ + -ip- . The 'ballpoint 
pen' model [19] is adopted wherein the spacetime metric describing the string's interior region p < p^ is 
matched smoothly at p = p^ to the spacetime metric describing the string's exterior region p > p^. The 
Tamm medium representing this cosmic spinning string is characterized by the constitutive relations [3] 



D{r, t) = eo7(a;, y) • E{r, t) - y/e^ r{x, y) x H{r, t) 
B{r, t) == pal{x, y) ' H{r, t) + ^^i^T^ r(a;, y) x ^(r, t) 



(1) 



in SI units. The scalar constants eo and po denote the permittivity and permeability of vacuum in the 
absence of a gravitational field, while the 3x3 dyadic 
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are expressed in terms of the real-valued constants A and a, which are related to the energy density and 

angular momentum respectively. As we discussed previously [3], the choice A = {ir/ps) ensures that the 

constitutive parameters of the Tamm medium coincide with those of vacuum (in flat spacetime) in the limits 

p — > and cxj. In keeping with our earlier study [3], we set a = 1 for our numerical studies of a spinning 

string, and a = for a non-spinning string. 

Consistently with the geometric-optics approximation, we consider quasi-planewave electric and magnetic 

fields of the form [20] 

E{L,t) = 'R-c{EJr)cicp[i{kok-r-ujt)]} ] 

}, (6) 

Kir, i) = Re { H„{r) cxp [i {k„k ■ r - cut)]} J 



where E_g{r) and iLoiL) a-re complex- valued amplitudes, a; is the angular frequency, and the vacuum 
wavenumber fco = ui,^eofi„. The relative wavevector k in the quasi-planewave representation (6) is a function 
of r but for later convenience we omit explicit reference to r. Combining eqs. (1) and (6) with the source-free 
Maxwell curl postulates yields 

[V (fc • r) - r(x,y)] X E,{r) - ,/^7(2:,2/) • ^o(l) = -^ V x E,{r), (7) 

[V(fc-r)-r(x,2/)]x^„(r) + ./^7(x,y)-^„(r) = -^Vx^„(r). (8) 

V Mo — *"'o 

In the geometric-optics regime, the constitutive parameters are assumed to vary very slowly over the distance 
of a wavelength. Accordingly, the right sides of eqs. (7) and (8) are set to zero, and V (A: • r) w fc. Thus, 
eqs. (7) and (8) simplify to [21] 

^^det2{x,y)-p-2ix,y)-p 1 + pp- gix^y)^ ■ E„{r) ^ 0, (9) 

where p = k — T_{x, y) and / denotes the identity 3x3 dyadic. The dyadic enclosed in braces on the left side 
of eq. (9) is required to be nonsingular, in order for nontrivial solutions to exist. This leads to the dispersion 
relation [21] 

n = det2{x,y)-p-2{x,y)-p = 0, (10) 

from which the magnitude k of the relative wavevector k may be deduced. In fact, the fc-roots of the 
dispersion relation (10) — corresponding to fc = fc {sm9 cost/), sin^ sini/), cos 6) at r = (pcostr, p sin cr, z) — 
may be expressed as 



-b ± Vb^ - 4:ac ,_,_,, 



where the coefficients 



a 



A^ cos^ 9 + [A^ cos^ (ct - 0) + P^ sin^ (ct - </>)] sin^ 9, (12) 

b ^ 2Mpsm9 sin{a~<j>), (13) 

c = M'^-A'^. (14) 

Although eq. (10) is a quadratic equation in fc, it yields only one independent root, in consonance with 
vacuum being unirefringent [22] . 

In general, k can be complex-valued with non-zero imaginary part. However, regimes where Im {k} ^ 
correspond to evanescent waves and arc therefore excluded from our study of ray trajectories. It is 
illuminating to characterize the phase space corresponding to Im {k} ^ 0. The discriminant term in eq. (11), 
namely, 

6^ - 4ac = 4 (m'^p^ sin^ 9 sin^ (a - cj)) 

- {A^ cos^ 9 + [A^ cos^ {<J -(!))+ P^ sin^ {a - <j>)] sin^ 9} [M^ - ^^) ) , (15) 

can only be negative- valued when A^ — Al"^ < 0. Therefore, the regime where evanescent waves arise coincides 
exactly with the regime where the spacetime of a cosmic spinning string can support closed timelike curves 
[3, 23]. This regime is illustrated in Fig. 1, wherein the directions of A: for which Im{A;} ^ are depicted at 
locations along the x axis. For reasons of symmetry, only the directions in one octant of the unit sphere need 
be shown. We sec that the Im{fc} 7^ regime is restricted to the range 0.1 < (p/ps) < 3. Within this range, 
Im{fc} 7^ for most k directions, but for k directed tangentially k remains real-valued. It is clear from eq. 
(15) that for a non-spinning string (i.e., M ~ 0), k is real-valued for all k directions, at all locations. 



3 Ray trajectories 

The scalar quantity % introduced in eq. (10) serves as a convenient Hamiltonian function for our ray-tracing 
study. The ray trajectories are parameterized in terms of r via rij) , and the relative wavevector is likewise 
parameterized as ^(t) . The coupled differential equations [24, 25] 



ar - 

f = -^^^ 
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(16) 



govern the trajectories of the light rays. Herein the shorthand Vi^ = {d/dvx, d/dvy, d/dvz) for v = 
{vx,Vy,Vz) is adopted. 

Before proceeding, we confirm that the direction of a ray trajectory, i.e., V k'H, coincides with the direction 
of energy flux. We know from earlier studies with a general Tanim medium that the time-averaged Poynting 
vector may be expressed as [26] 

(P)t-/32(a;,y)-p, (17) 

where the scalar /3 is positive-valued provided that "f{x,y) is either positive- or negative-definite. From 
eq. (2), the eigenvalues of ^{x,y) are {A/ p,A/ p,p/A}] hence /3 > 0. By direct vector differentiation of the 
definition of "H provided in eq. (10), we find V^H = 2-f{x,y) ■ p . Therefore, the ray trajectories extracted 

from eqs. (16) are indeed parallel to the direction of energy fiux. 

Upon the specification of appropriate initial conditions r(0) and fc(0), the system (16) can be solved for 
r^r) — and k(j-) — by standard numerical methods, such as the Runge-Kutta method [25]. For a spinning 
string, two examples of ray trajectories are plotted in Fig. 2. In the first example, the rays start at locations 
in the xy plane, with k(0) directed along (—1, 0, 0). In this case the ray trajectories remain in the xy plane. 
In the second example, fc(0) is directed along (—1,0,-1) and the ray trajectories are not restricted to one 
plane. In both examples, the ray trajectories skirt around the string's boundary, never crossing it. 

Let us now consider further whether it is possible for ray trajectories to cross the string's boundary. In 
Fig. 3, trajectories are presented for rays which start close to the string boundary, at r(0) = (1.1, 0, 0) ps and 
(0.9, 0, 0) Ps- Regardless of the direction of fc(0), we find that rays which start outside the string's boundary 
remain outside, and those which start inside the string's boundary remain inside. We therefore conclude 
that rays cannot cross the string's boundary, in either direction. Note that a ray starting right at the string's 
boundary cannot be considered, as the geometric-optics approximation is not valid at this location. 

4 Invisibility 

The ray trajectories presented in Fig. 2 are reminiscent of those associated with metamaterials which are 
currently being investigated for cloaking applications [18, 27, 28, 29]. That is. the region inside the spinning 
string's boundary is not visible to a distant observer. However, this self-cloaking cloaking effect is not perfect 
as the ray trajectories reaching a distant observer are distorted somewhat by the string. 

In order to explore this matter further, let us consider a non-spinning string. The ray trajectories for a 
non-spinning string — corresponding to those presented in Fig. 2 for a spinning string — are presented in 
Fig. 4. We see that in a radial plane, the non-spinning string acts as an excellent cloak for itself: the interior 
of the string is hidden to a distant observer and there is minimal distortion to the ray trajectories reaching 
an observer more than a few radiuses away from the string. For rays traversing the string's neighborhood 
at an oblique angle to the z axis, there is a small degree of ray distortion close to the string — but this 
distortion diminishes as the distance from the string increases. 



5 Closing remarks 

The Tamm medium provides a convenient setting for simulating the passage of hght through the spacetime 
associated with a cosmic spinning string. Based on a geometric-optics study, we have found that 

• ray trajectories do not cross the string's boundary; i.e., rays which start outside the string remain 
outside, whereas rays which start inside the string's boundary remain trapped there; 

• the regions of the spinning string's spacetime which support cfosed timehke curves correspond to regions 
in the phase space of the Tamm medium which support evanescent waves; and 

• a spinning string acts as an imperfect cloak for itself, while a non-spinning string cloaks itself almost 
perfectly. 

These findings may well have far-reaching astronomical consequences, especially for those attempting to 
observe cosmic strings optically. 
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Figure 1: Maps illustrating the directions oi k = k ( k^, ky, kz 1 for which Im {fc} = (blue) and Im {fc} ^ 
(red), at locations along the x axis with x e {0.1, 0.2, 0.9, 1.1, 2.7, 3.0} p^- 





Figure 2: 2D and 3D examples of ray trajectories for a spinning string. Top: Rays start at equally spaced 
locations along the line a;(0) = (40, :/, 0) ps with —3 <v < —1, and fc(0) directed along (—1,0,0). Bottom: 
Rays start at equally spaced locations along the line x(0) = (40,i^, 39)ps with —3.9 <v < —1.9, and fc(0) 
directed along (—1,0,-1). 
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Figure 3: Examples of 2D ray trajectories that start at locations close to the spinning string's boundary: 
r(0) = (1.1,0, 0) ps (top) and (0.9, 0, 0) ps (bottom). In both cases fc(0) is oriented at equally spaced angular 
directions in the xy plane. Note that for both locations approximately 70% of the possible orientations of k 
correspond to evanescent waves and are therefore not represented here. 
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Figure 4: As Fig. 2 except that the string is not spinning and — 1 < i^ < 1. 
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